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Abstract. Assuming the magnetic hypothesis for the for- 
mation of warps, we deduce a very simple formula for the 
warp curve in an idealized scenario. According to this for- 
mula the warp rises as the third power of radius in the 
innermost warped region, reaches the maximum slope at 
intermediate radii and has an asymptotic slope coincident 
with the direction of the extragalactic magnetic field. In 
most cases, however, the galaxy's limited size prevents the 
observation of the full curve. Even though the model is 
highly simplified, it basically reproduces real warp curves, 
in particular the 21 cm warp curve of NGC 5907. If the 
magnetic model were considered to be correct, the fit- 
ting of warp curves could allow rough estimations of the 
strength and direction of the magnetic field. We also pro- 
pose a magnetic field distribution for the outermost part 
of the galaxy, matching the boundary conditions of be- 
ing azimuthal inside, and constant at infinity. We use this 
magnetic field distribution to show that the assumptions 
made to obtain the warp curve with our simple model 
cannot introduce important errors. 
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1. Introduction 

The possibility of warps being induced by extragalactic 
magnetic fields (Battaner et al., 1990) remains a tempt- 
ing explanation of this common feature of most spiral 
galaxies. Without intending to enter a discussion about 
its validity against other current models, some observa- 
tional facts may produce a renewed interest in the mag- 
netic model: a) Binney (1991) and Combes (1994) argued 
that the required magnetic fields were too high for the in- 
tergalactic medium. However, Kronberg (1994) reviewed 
evidence for magnetic fields of about 1-3/xG to be com- 
mon in intra-cluster media. Later, Feretti et al. (1995) 
have found B < 8.3/xG in the Coma cluster, about one 
order of magnitude higher than the equipartition value. 
This is a critical point, as magnetically driven warps re- 
quire magnetic strengths of this order of magnitude. Fur- 



ther confirmation of high intergalactic fields would sup- 
port our hypothesis. We are not aware of any published 
alternative interpretation of the available observations of 
intergalactic magnetic fields, b) One of the most interest- 
ing possibilities is that of warps being normal modes of 
oscillation of the galactic disk (Sparke and Cassertano, 
1988). But when the back reaction of the halo is taken 
into account this model seems to fail (Dubinski and Kui- 
jken, 1995; Kuijken, 1997). c) Zurita and Battaner (1997) 
have shown the tight alignment of the warps of the three 
largest spirals in the Local Group. Such an alignment is 
easily understood in the magnetic model, and probably 
also in the model by Kahn and Woltjer (1959). 

In this short paper we obtain a simple formula describ- 
ing the warp curve [mean z versus x] when the magnetic 
hypothesis is adopted. The coordinate x is contained in 
the inner unwarped disc plane and is perpendicular to the 
line of nodes, assumed to be untwisted for simplicity. Very 
simple models for idealized scenarios are specially appro- 
priate for exploratory hypotheses and remain useful even 
when more sophisticated computations are available. 

2. The shape of magnetically induced warps 

According to Battaner et al. (1990), the vertical magnetic 
force which induces the warp is of the order ByB z /8irL, 
where B y is the component of the extragalactic field in 
the equatorial plane of the galaxy perpendicular to the 
line of nodes, B z is the component in the direction of the 
rotation axis of the galaxy and L is an undetermined char- 
acteristic length. To deduce this expression, we assume a 
constant mean value of the magnetic field strength in a 
large characteristic length connecting the inner disc and 
the intergalactic medium. This seems to be an extremely 
simple assumption, but it will be shown later that it does 
not introduce a quantitatively important error. Let us as- 
sume in this simplified model that these quantities B y , 
B z , and L are constants in the region of interest. 

The characteristic length L would depend on the de- 
gree of ionization which in turn depends on the galacto- 
centric radius. Let us however assume that at any radius 
the degree of ionization is large enough as to assure in- 
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finite conductivity and frozen-in magnetic field lines. In 
the inner part of the disc this force is negligible compared 
with the gravitational force, but it becomes increasingly 
important towards the outer parts of the disc. 

The expression is valid for a point in the equatorial 
plane of the galaxy. For a point above this plane the ver- 
tical force becomes 



F z = 



B 2 sin(2(q-/?)) 
WttL 



cos/3 



(1) 



where tan [3 — z/x and tan a is the slope of the direction 
of the magnetic field in a [a;, z] diagram. It is easily derived 
that this expression is equivalent to 



B 2 x{x 2 - z 2 ) . , n , 



2zx 2 



cos(2a) 



(2) 



where r = (x 2 + z 2 ) 1 / 2 is the radial coordinate in the [a;, z] 
plane 

In this simple model we assume that the gravitational 
potential is that of a point mass in the centre of the galaxy, 
as a simplifying assumption for the outer part of a disk not 
embedded in a massive halo. This potential has already 
been used before, for instance by Cuddeford and Binney 
(1993). Equilibrium in the vertical direction gives 

dp 



dz 



pGMr^z 



B 2 



WttL 



x(x 2 ~ z 2 ) 



sin(2a) 



2zx 2 



cos(2a) (3) 



A full solution for the distribution of the gas in the 
combined magnetic-gravity force field would reproduce the 
whole geometry of the warp. This full solution is beyond 
the scope of this paper. In any case, it would be interest- 
ing to, simply bu precisely define the warp curve. In the 
inner unwarped region, ^ = at z = 0, in the galactic 
plane. Let us therefore define the warp curve as the locus 
of points where dp/dz — 0. We adopt an exponential law 
for the disc density with length scale R, i.e. p — poe~ x ' R . 
Using R as length unit for x and z, we obtain for the warp 
curve 



[-e~ x - 2kx 2 cos(2a) + {(e~ x + 2kx 2 cos(2a)) 2 
+4fcV sin(2a)} 1/2 ] (2kx sin(2a))" 1 



where 



k = 



R 2 B 2 



(4) 



(5) 



16irLp GM 

is one of the adjustable parameters which compares the 
extragalactic magnetic energy density with the gravita- 
tional energy density. The other adjustable parameter is 
a which specifies the direction of the extragalactic mag- 
netic field. The warp curve is therefore defined with just 
two free parameters: k and a. 

For small values of £ a series expansion gives 



which is a very simple expression for small warps. This 
simple formula illustrates the fact that the maximum ef- 
ficiency in producing warps is obtained for a = 45°, as in 
Battaner et al. (1990) 

For very large values of x, we obtain 



z = tan(a)x 



(7) 



k sin(2a)x 3 



(6) 



i.e. the slope of the warp curve matches the direction of 
the magnetic field at large radii. 

In Fig. [j] we plot the obtained curves for k = 
l(T 3 ,l(r 4 ,Kr 5 ,lfr 6 ,lCT 7 for a = tt/4 and a = tt/6. 
The curves for a — 7r/4 seem to be unrealistic at first 
glance. Note, however, that in practice values at x » 6 
are unobservable (or the galaxy simply does not exist at 
these radii). For instance if the [a — 7r/4, k = 10~ 5 ] curve 
is truncated at x ~ 6, the obtained warp curve becomes 
quite familiar. We reproduce the curve for larger values of 
x and z in order to see the region where the slope becomes 
equal to the direction of the extragalactic magnetic field, 
which probably takes place for radii far from observational 
capacities or galaxy limits. 

It is worth noting that, for a < tt/4, there is a change 
in slope. It is higher at intermediate regions before reach- 
ing its asymptotic value tana. This is a common feature 
of real warps, and can even be directly observed in the 
early contour maps of NGC 5907 and NGC 4565 by San- 
cisi (1976). 

Figure |^ reproduces the observational curve for one of 
the best known prominent warps, in NGC 5907, adopted 
from Sancisi (1976). We cannot exclude that the warp of 
this galaxy is due to other mechanisms, but we choose 
this warp because it is one of the most representative and 
is studied very often. This figure also reproduces a fit- 
ting to the model, with parameters k = 8 x 10~ 5 and 
a = 20°. It is not straightforward to deduce the mag- 
netic field strength from the value of k, mainly because L 
is an equivalent undetermined quantity. For L = 1 kpc, 
pa = 1.7 x 10~ 24 g cm" 3 , M = 2 x 10 11 M© and 
R = 4 kpc, it is obtained that values of k in the range 
10 -4 to 10 -6 correspond to field strengths between 3 pG 
and 0.3 pG, in agreement with reported measurements by 
Kronberg (1994). Probably better results would probably 
be obtained with more realistic calculations, but given the 
present, still exploratory, character of the magnetic model 
it is preferable to deal with idealized systems. The notice- 
able fitting of the NGC 5907 warp curve suggest that the 
magnetic model and the approximations considered here 
are not unreasonable 



3. Magnetic field distribution in the outer disc 

We have studied some magnetic effects in the density dis- 
tribution of the outermost disc without considering in de- 
tail the magnetic field distribution in this region. This was 
made on one hand because this distribution is unknown, 




NGC 5907 




Fig. 2. Warp curve for NGC 5907. Open symbols are 
experimental data adopted from Sancisi (1976) and the 
solid line corresponds to the model for k = 8 x 10~ 5 and 
a = 20°. 



and on the other hand to produce an exploratory simpli- 
fied model. Now, we propose a magnetic field distribution 
and use it to estimate the error of our assumption for the 
specific case of a = 45°, small warps and large galacto- 
centric radii. 

As no radio continuum measurements are available to 
deduce the magnetic field in this region, we must just as- 
sume it. For this task we have not so many degrees of 
freedom, because our field must fulfill two restrictions: a) 
Boundary conditions: it should be azimuthal (B^) at the 
inner radius L\, and should become constant at infinity, 
b) The field must also have vanishing divergence. We thus 
solve the general problem of connecting both, the disc 
magnetic field and the extragalactic magnetic field in a 
way that could be adopted in other problems not associ- 



ated with the particular case of warps. Outside the galaxy 
the field could have any direction. To consider a general 
case we take B x = B z = B^j \/2 = A outside, but B y = 0, 
due to a proper choice of the y axis. Let us define: 



u 



A * 2 +f . 2 2 , 

—e ^^(y 2 - a 1 ) +A 



B2x = B u 



B> 



2y 



- & 



- 1/2 L 2 -(x 2 + y 2 )^ 2 



-1/2 



Li 



(x 2 + y 2 ) 1 ' 2 



Z/2 — L\ 



(8) 
(9) 

(10) 
(11) 



We then consider three regions: 



1. Inner disc, for < r < L\. We do not propose any field 
distribution for this region, as our model only considers 
the periphery of the galaxy. 

2. Outer disc, for L\ < r < L2, taking L2 as the radius 
where the galaxy's edge is assumed to be. In this re- 
gion, 



B, f - 



B\ x + B 2x 



£>„ — B 



B- 



2y 



(12) 
(13) 



3. Nearby intergalactic space, for r > L 2 . We propose in 
this region, 



Br,, — B 



B y - B ly 



(14) 
(15) 
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Fig. 3. Magnetic field distribution of the galactic disc in 
the XY plane for B v = 3/iG, A = 0.2(iG, L 1 = 2,L 2 = 'l 
and a = 2 (taking the radial scale length as unit). 



Fig. 5. Magnetic field lines in the outer disc and in the 
intergalactic space close to a spiral galaxy. The warp plane 
is the y = plane. 



Plane Y=0 




Fig. 4. Magnetic field distribution of the galactic disc in 
the Y=0 plane for B v = 3/zG, A = 0.2^G, L\ = 1, L 2 = 3 
and <j = l. 



It can be easily checked that the divergence of this two- 
dimensional field is zero. In Fig. |3|, we plot this magnetic 
field distribution. The magnetic field strength along x is 
plotted in Fig. |]. This distribution has no jump (only the 
first derivative is not a continuous function). In Fig. |] 
there is an apparent jump at L\, but it only means that 
we are ot interested in the inner disc. 

Some of the formulae above would become clearer in 
polar coordinates, but cartesian coordinates are the nat- 
ural ones outside the galaxy, and a choice covering the 
whole region is preferable. 



We must complete our choice by defining the z compo- 
nent. We assume that B z is independent of z, in order to 
conserve the vanishing divergence. Then, for L\ < r < L 2 

B ,= v^ + / )1, ;- L ' <i6) 

For r > L 2 we adopt B z = A. 

In this case, in the direction of the warp plane (per- 
pendicular to the line of nodes and to the galactic plane) 



dB z A 
^ {y = ° )= L- 2 



(17) 



as in Battaner, Florido and Sanchez-Saavedra (1990). The 
three-dimensional picture of the magnetic field distribu- 
tion can be appreciated by combining figure @ and eq. 
(|l6|). The result is depicted in figure (^), where the inner 
azimuthal field lines have not been plotted. 

Let us now consider a second model, with the simplify- 
ing assumption of a = 45° and z warp <ti x warp , but taking 
into account the above proposed magnetic field configura- 
tion. 

The vertical force will be: 



F, 



1 dB z 
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47T OX 
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(18) 
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And for y=0 (warp plane) 



4irL? 



The equation of vertical equilibrium becomes 



^ + P GMr- 3 z = 
oz 



B n 



8ttL 2 



1 — e 20 



(19) 



(20) 



To calculate the warp curve, we use || = as before. We 
now find for z small (so that r ss x) 



z = 



1 



(21) 



8ttL 2 GMpo 

With the same assumptions the first model would give 

= Boo « , . 

8irL 2 GMp Q V ; 

We note, however, that the last exponential function in eq. 

( |2l| ) is not very important. For x — a, we have — 
0.61; for x = 2a, we have 0.14; for x — 3a, we have only 
0.01. Therefore, we conclude that the use of the first model 
is fully justified. 
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4. Conclusions 

In order to study the influence of intergalactic magnetic 
fields on the warping of galactic discs, it is necessary to 
know their distribution in the outermost disc, B(x,y,z). 
At present, no observations are available to determine the 
field at very large galactocentric radii. However, we have 
proposed a field distribution which should not differ very 
much from the real one, as the boundary conditions and 
the condition V ■ B are in practice very restrictive. Mag- 
netic field lines must be very similar to those presented in 
our figures [|, || and ||. 

We have shown that the detailed knowledge of the 
magnetic field lines in the periphery of the galaxy is rel- 
atively unnecessary to predict reasonable warp curves. 
More sophisticated calculations should be required in the 
future, but considering the present state of our under- 
standing of warps, a simplified model has more advan- 
tages. The magnetic model of warps provides curves in 
good agreement with real warps. Under the magnetic hy- 
pothesis, the fitting of the warp geometry could provide 
information about two important cosmological parame- 
ters: the direction and the strength of the large scale ex- 
tragalactic magnetic field. 
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